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Abstract: Pholosystem 11 (PSII) continuously allracls Ihe allenlion of researchers aiming to unravel the riddle of its func- 
tioning and efficiency fundamental for all life on Earth. Besides, an increasing number of biotechnological applications 
have been envisaged exploiting and mimicking the unique properties of this macromolecular pigment-protein complex. 
The PSII organization and working principles have inspired the design of electrochemical water splitting schemes and 
charge separating triads in energy storage systems as well as biochips and sensors for environmental, agricultural and in- 
dustrial screening of toxic compounds. An intriguing opportunity is the development of sensor devices, exploiting native 
or manipulated PSII complexes or ad hoc synthesized polypeptides mimicking the PSII reaction centre proteins as bio- 
sensing elements. This review offers a concise overview of the recent improvements in the understanding of structure and 
function of PSII donor side, with focus on the interactions of the plastoquinone cofactors with the surrounding environ- 
ment and operational features. Furthermore, studies focused on photosynthetic proteins structure/function/dynamics and 
computational analyses aimed at rational design of high-quahty bio-recognition elements in biosensor devices are dis- 
cussed. 

Keywords: Molecular docking, molecular dynamics simulations, Photosystem II, plastoquinone, plastoquinone binding site, 
protein dynamics. 



1. STRUCTURAL OVERVIEW OF PSH CORE COM- 
PLEX 

Photosynthetic organisms have the unique ability of en- 
trapping and using solar energy to oxidize water to molecular 
oxygen through the highly structured multi-subunit protein- 
chlorophyll complexes of photosystem I (PSI) and photosys- 
tem 11 (PSII), arranged in supra-molecular assemblies em- 
bedded in thylakoid membranes. These structured assembhes 
are organized in a way to ensure supreme energy conversion 
efficiency and effective diffusion of small molecules (plas- 
toquinone) interconnecting the two photosystems functions 
in a fine-tuned electron transport chain and leading to the 
synthesis of ATP and reducing equivalents. PSII stands out 
for its remarkable ability to act as a water:plastoquinone oxi- 
doreductase, extracting electrons from water in physiological 
conditions. The outstanding features of PSII have motivated 
numerous research efforts aimed to clarify its structural and 
functional relationships, which until recently have been hin- 
dered by the lack of a high resolution X-ray crystal structure. 
Our current understanding of the PSII struc- 
ture/function/dynamics rationale has been largely acquired 
by comparative studies and assumptions with the simpler and 
better-known purple bacterial reaction centre (bRC) 
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[for review on bRC see 1, 2]. In the last decade, the crystal 
structure resolution at 2.9-3.5 A has exponentially improved 
the knowledge of atomic details of dimeric and monomeric 
forms of PSII core complex (PSIIcc) of the thermophilic cya- 
nobacterium Thermosynechococcus elongatus [3-7]. More 
recently, Umena and coworkers provided the highly resolved 
1 .9 A structure from Thermosynechococcus vulcanus [8] . This 
structure features a rich collection of atomic details among 
which a cofactors arrangement made up by 35 chlorophylls, 2 
pheophytins, 11 (3-carotenes, 2 plastoquinones, 2 heme irons, 
1 non-heme iron, 4 Mn atoms, 3-4 Ca atoms (among which, 
1 is in the Mn4Ca-cluster), 3 CI" ions (among which, 2 are in 
the vicinity of the Mn4Ca-cluster), and 1 (bi)carbonate ion. 
From the electron density map it was possible for the first 
time to locate all the atoms building up the Mn4Ca05-cluster 
and to define the precise position of all the metal ligands emd 
of the substrate wata" molecules. In this high-resolution struc- 
ture the primary plastoquinones (PQs) Qa and Qb occupy po- 
sitions very similar to those observed in previous crystallo- 
graphic studies. However, the third plastoquinone molecule 
observed in the structure solved by Guskov and coworkers [5], 
the so called Qc, has not been observed in this last structure 
casting some doubts on the real physiological function of Qc. 
These structural data together with the numerous biochemical 
and spectroscopic methods, and computational analyses con- 
tribute to solve the enigma of PSH structure/function/ 
dynamics specific correlation. 
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PSIIcc consists about 20 different polypeptide chains, 
with the majority in the form of trans-membrane and a- 
helical subunits. Among these, PsbA (Dl) and PsbD (D2) 
subunits, with five trans-membrane hehces (TMH) each, 
form the architecture of PSII RC and host the major redox- 
active cofactors organized in two quasi symmetrical 
branches: two main chlorophyll a (Chla), PDl and PD2, two 
accessory Chla, ChlDl and ChlD2, two peripheric Chla, 
ChlzDl and ChlzD2, two Pheophytin a, PheoDl and 
PheoD2, and two PQs, Qa and Qb (Fig. 1). On the stromal 
(cytosol in prokaryotes) side of the complex between the 
primary and secondary PQs, a non-heme iron interacting 
with a bicarbonate anion has been also revealed. It is electro- 
statically associated to the two PQs, optimizing charge dis- 
tribution and H-bonding of the QA-Fe-Qe bridge and con- 
tributing to the stability of PQs when reduced to semiqui- 
nones (SQ) [9, 10]. Two antenna proteins, PsbB (CP47) and 
PsbC (CP43), are located in close proximity of the D1/D2 
heterodimer. Their structure consists of six TMHs each and 
large extrinsic loops at the lumenal side. They host the ma- 
jority of PSIIcc chlorophyll molecules (16 and 13 Chla, re- 
spectively) and, besides light harvesting, are involved in the 
assembly and stabilization of the Oxygen Evolving Complex 
(OEC) active conformation [8, 11]. Although the resolution 
of the PSII crystal structure at 1.9 A provided an improved 
picture of the overall configuration of the OEC Mn4Ca05 
cluster, and atomic distances of the Mn-Mn and Mn-Ca in- 
teraction network, the structural changes occurring during 
catalysis are still under investigation. Different X-ray ab- 
sorption spectroscopy techniques and quantum mechani- 
cal/molecular mechanics calculations are providing further 
insights into the relative position of one of the bridging oxy- 
gen in the Mn4Ca05 cluster, that are crucial for a detailed 
understanding of the 0-0 bond formation mechanism during 
the water oxidation reaction [12, 13]. An overview of the 
kinetics and thermodynamics of water oxidation and new 
knowledge on biogenesis and assembly of OEC has been 
critically addressed recently by Vinyard and coworkers [14]. 
Among the major intrinsic proteins necessary for PSII func- 
tion are the TMHs of PsbE and PsbF, corresponding to a- 
and ^-subunit of the heme-containing protein cytochrome 
^559 (Cyt.bssg), which most probably takes an active part in 
PSII photoprotection mechanisms. However, its involvement 
in PSII forward electron transport seems unlikely, despite its 
location close to the secondary quinone binding site [15, 16]. 

Important molecular players of PSIIcc are also carotenoid 
(Car) molecules; the position of 12 p-carotenes per PSIIcc 
monomer has been assigned in the PSII structure refined at 
2.9 A resolution [5]. It has been proposed that these poly- 
isoprenoid compounds can play crucial roles in light- 
harvesting, excitation energy transfer regulation, Chi triplet 
states quenching, singlet oxygen scavenging, and structure 
stabilization [17]. Crystallographic studies have also re- 
vealed up to 25 integral lipids associated to each PSIIcc 
monomer, including 11 monogalactosyldiacylglycerol 
(MGDG), 7 digalactosyldiacylglycerol (DGDG), 5 sulfoqui- 
novosyl-diacylglycerol (SQDG), and 2 phosphatidylglycerol 
(PG), having almost full compositional and spatial overlap in 
PSII monomeric and dimeric form [5, 7, 8]. PSII lipid con- 
tent and composition denote the typical thylakoid membrane 
distribution and are characterized by spatial heterogeneity: 



the head group has a stromal (cytoplasmic in prokaryote) 
orientation for the negatively charged PG and SQDG, lume- 
nal for DGDG and bilateral for MGDG. PSII integral lipids 
behave more as a dynamic multifunctional cofactor into 
overall functioning of the complex than as a simply hydro- 
phobic surrounding matrix. Lipid molecules, in fact, actively 
interact with the Mn-cluster components and stabilize its 
assembly; contribute to a different extent to define the con- 
figuration of approximately 1/3 of Chla and 3/4 of Car mole- 
cules in PSII, and last but not least, largely contribute to the 
PSII monomer-monomer interactions and complex dimeriza- 
tion [18, 19]. Integral lipids of PSIIcc play a special role also 
in defining the large internal cavity for the PQ/PQH2 ex- 
change, representing the ideal flexible scaffold to host the 
hydrophobic PQ isoprenoid residues, including the recently 
identified PQ molecule, Qc. The cavity has two openings 
directed towards the membrane heart, approximately in the 
centre of the membrane thickness, defined as channel I and 
channel II [5]. Channel I is bigger than channel II, lodges 
Qc, and opens between TMHs of PsbJ protein and Cyt.bssg, 
while the recently described channel II accommodates Qb 
phy tyl tail and is formed by TMH-d and e of D 1 , TMH-a of 
D2 and TMH of PsbF (^-subunit of Cyt.bssg)- 




Fig. (1). Schematic representation of PSIIcc main protein compo- 
nents and electron transfer cofactors. The reaction centre (D1/D2) 
heterodimer is reported in light (D2) and dark (Dl) yellow, respec- 
tively; core antenna proteins CP47 and CP43 are in green; the Oxy- 
gen Evolving Complex (OEC) - in blue. PSIIcc associated cofac- 
tors: the primary (PD1/PD2) and accessory (ChlDl/ChlD2) Chla 
molecules are in cyan, pheophytin a (PheoDl/PheoD2) are in gold 
brown, plastoquinones (Qa, Qb and Qc) are in yellow, the heme of 
cytochrome b^sg (Cyt. b^sg) is colored green and the non-heme iron 
in grey. 

2. PSII PLASTOQUINONES 

The two chemically identical PQ molecules, hosted into 
the Qa and Qb binding sites of the RC heterodimer, have the 
distinctive role of stabilization of primary charge separation 
and conduction of electrons towards the photosynthetic elec- 
tron transport chain, respectively. The existence of a third 
plastoquinone, Qc, in the PSII core complex has been sup- 
ported by Cyt.Z)559 photoreduction and redox titration ex- 
periments and lately identified by the 2.9 A resolution crys- 
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tallographic study [5, 20-22]. The role of this quinone mole- 
cule in the PSII function is still under debate and will be dis- 
cussed below. The binding energies for Qa, Qb and Qc ob- 
tained by fragment molecular orbital (FMO) calculations 
based on the 2.9 A X-ray structure are -56.1 kcal/mol, -37.9 
kcal/mol and -30.1 kcal/mol, respectively [23]. The FMO 
calculations indicated a decrease of the relative contribution 
of the PQ/protein matrix interactions to the total PQ interac- 
tion energy in the order Qa>Qb>Qc with concomitant in- 
crease of the relative contribution of surrounding lipids and 
cofactors (Table 1). Thus, the FMO results drew a correla- 
tion between the increasing of PQ mobility and the growing 
relative contribution of lipids and cofactors interaction. 

2.1. Plastoquinone at Qa Binding Site 

The plastoquinone at Qa site tightly interacts with the RC 
protein matrix, does not undergo protonation, and is not ex- 
changeable. It functions as a single electron recipient and 
mediates the double reduction of the second quinone mole- 
cule at the Qb binding pocket. In contrast, the double re- 
duced and protonated Qb moves away from its binding site 
and is promptly replaced by a new PQ molecule (as well as 
by artificial quinones or herbicides), leading to consider this 
molecule more as a substrate rather than a tightly bound co- 
factor. Although chemically identical, Qa and Qb demon- 
strate distinguishable redox potentials, necessary for the po- 
tentiation of PSII forward electron transfer and the fine tun- 
ing of the backward recombination reactions [24, 25]. The 
different functional and redox-properties of the PSII plas- 
toquinone molecules are largely determined by the interac- 
tions with the surrounding environment (H-bonding, hydro- 
phobic interfaces, 7t-stacking), and the conformation of the 
individual binding sites and phytyl chains [4-6]. FMO quan- 
tum chemical calculations based on the 2.9 A X-ray structure 



of T. elongatus pointed out that the Qa phytyl tail is respon- 
sible for more than 60% of the total molecule interactions 
and Qa stabiUzation is mainly due to interactions with amino 
acid residues [23], Table 1 and Table 2. The position of the 
isoprenoid tail is restrained mostly by van der Waals contacts 
with TMH-^i of D2, TMH-a of Dl, PsbL and PsbT and the 
accessory ChlDl (Fig. 2). The 2.9 A resolution structure 
revealed an additional stabilization of the Qa molecule due 
to the interaction of the terminal unit of the isoprenoid tail 
with 7t-electron clouds of the phenyl rings of Dl-Phe52 and 
PsbT-PhelO [5]. Qa head is fixed by two H-bonds, formed 
by the keto-oxygen atoms with the D2-Phe261 backbone 
amide group and the side chain N of D2-His214, and by n- 
stacking with D2-Trp253 indole ring, occurring in an offset- 
stacked manner [4, 5]. The static picture of Qa interactions 
within its binding niche obtained by X-ray analyses has been 
further complemented by spectroscopic studies. Simultane- 
ous 'H and "^N Two-Dimensional Hyperfine Sublevel Corre- 
lation spectroscopy shed new light on the strength and orien- 
tation asymmetry of H-bonding of Qa anion radical, estab- 
lishing that when Qa is in its reduced form, the H-bond in- 
volving D2-Phe261 amide group is significantly stronger 
than the one attributed to D2-His214 [26, 27]. Similarly, 
computational studies provided stronger interaction energies 
of Qa with D2-Phe261, D2-Thr217 and D2-Trp253 than with 
QA-D2-His214 [23], (Fig. 2, left bottom panel). 

Important properties of the primary PQ molecule, proven 

in bacteria, cyanobacteria and higher plants, is the existence 
of the one-electron redox couple Qa/Qa' in two conforma- 
tions with highly different redox potential [28]. The well- 
recognized ability of herbicides binding at Qb site to increase 
or decrease the susceptibility of PSII to photodamage was 
related to herbicide-induced changes at midpoint potential of 
the Qa [29-31]. Recent studies have strengthened a previ- 



Table 1. Distribution of interaction energy of the PSII plastoquinone within Qa, Qb and Qc binding niche, according to Hasagawa 
and Noguchi [23]. The relative contribution of the protein matrix, and lipids and cofactors, respectively, are calculated as 
a percentage of the total interaction energy of the corresponding PQ molecule. The head/tail distribution of the total PQ 
interaction energy, as well as with the protein matrix and lipids and cofactors are calculated in a similar manner. The 
relative contribution of the PQs phytyl tail and head group for the total molecule interactions are also presented. Correla- 
tion between the increased PQ mobility (Qa>Qb>Qc) and reduction in the interaction with surrounding protein matrix is 
underlined. 



Factors Contributing to PQs Interactions 


Interaction Energy Distribution, % 


PQ at Qa Binding Niche 


PQ at Qb Binding Niche 


PQ at Qc Binding Niche 


Protein matrix PQs interactions 


100.0 


82.0 


56.0 


PQs head contribution 


38.0 


44.9 


8.8 


PQs tail contribution 


62.0 


56.1 


91.2 


Lipids and cofactors PQs interactions 


None 


18.0 


44.0 


PQs head contribution 




0.0" 


10.0 


PQs tail contribution 




100.0 


90.0 


Relative head contribution to total PQs interaction energy 


38.0 


35.6 


44.6 


Relative tail contribution to total PQs interaction energy 


62.0 


64.4 


55.4 



"The repulsive Qb interactions with Mpids and cofactors counteract to the attractive one. 
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ously found correlation [32, 33] between the Qa/Qa* niid- 
point potential changes and PSII donor side integrity level 
[28, 34]. The Qa low potential form is present in functional 
PSII reaction centres while donor-side Ca^"^ extraction or Mn 
cluster disassembly is associated with an increase (of approx. 
150 mV) in the Qa redox potential. Muh and coworkers [35] 
and Cardona and coworkers [36] have recently addressed in 
comprehensive reviews the physiological significance of the 
Qa redox shift in guiding a "safe" PSII charge recombination 
reactions and protecting PSII from photodamage during the 
assembly of the RC. In addition, an up-shift of the Qa mid- 
point potential by ca. h-38 mV was observed in T. elongatus 
in photo- oxidative stress conditions, when in the PSII RC the 
PsbA3 (Dl) isoform preferentially replaces the PsbAl (Dl) 
protein, expressed in optimal conditions. This finding evi- 
denced the capability of the Dl amino acid composition to 
influence the Qa redox potential, most probably affecting the 
QA-His214(D2)-Fe-His215(Dl)-QB interaction chain [37]. 
The rationale of the Qa redox potential switching in photo- 
protection should be viewed in light of the very recent study 
conducted by Boussac and coworkers [38], demonstrating 
faster recombination reaction in PsbA3 than in PsbAl PSII- 
RC, which could thwart 'O2 generation initiated by triplet 
state of the RC primary chlorophyll molecules. 

2.2. Plastoquinone at Qb Binding Site 

The site hosting the secondary quinone, Qb, is situated in 
a symmetry-related manner to Qa binding spot around the 
non-heme iron (Fig. 2, top panel). Because of its substrate- 
like behaviour it was not always possible to assign the elec- 
tron density found in the presumed Qb binding site to a PQ 
molecule [6]. The problem of the Qb binding site electron 
density inhomogeneity was surpassed by crystallization of 
the PSII core complex paired with a herbicide molecule [7] . 
However, in the 2.9 A structure it was possible to model 
seven out of the nine Qb isoprenoid units [5, 11]. In contrast 
to the secondary quinone acceptor in bRC, found in two dif- 
ferential positions called "proximal" £md "distal" [39], PSII 
Qb was located in only one position, corresponding to the 
"proximal" one in bRC. The Qb binding site is formed 
mainly by the Dl RC protein residues (Table 2). The keto- 
oxygen atoms of the Qb head group form H-bonds with Dl- 
Ser264 hydroxyl group and the side chain N of Dl-His215. 
In contrast to Qa, Qb has a possibility to form a third H-bond 
with the backbone amide group of Dl-Phe265 (2.84 A) [5, 
8], (Fig. 2, centre panel). Most probably, the latter is in- 
volved in the stabilization of the SQ state of the Qb molecule 
during the reduction/protonation process [36]. The Qb iso- 
prenoid tail, located into the recently described PQ/PQH2 
exchange channel II, is inchned in relation to the one of Qa 
[5]. Similarly to Qa, the contribution of the phytyl tail to the 
total Qb interaction energy is calculated to be more than 
60%, but differently from the Qa the surrounding protein 
matrix contributes for a smaller fraction of the binding en- 
ergy [23], Table 1. Based on the interaction strength of the 
Qa and Qb with the respective amino acid residues obtained 
by FMO calculations, Hasegawa and Noguchi [23] have 
questioned the role of the H-bonds established by the PQ 
keto-oxygen groups in the Qa and Qb binding. The authors 
observed more substantial contribution of the non-covalent 
interactions of Qa head group with D2-Phe261 and D2- 



Trp253 and Qb head group with Dl-Phc255 and Dl-Phe265 
to the total PQ binding energies than the above mentioned H- 
bonding. Instead, it has been unambiguously demonstrated 
the importance of the H-bonding patterns of PSII PQs in 
determining their redox potential [40, 41]. The relative con- 
tribution of PSII integral lipids and cofactors to the total Qb 
interaction energy is less than 20% and is mainly limited to 
the PQ tail (Table 1 and Table 2). The Qb phytyl chain is 
accommodated in channel II of the PSII internal cavity, con- 
necting the Qb binding site with the membrane PQ pool, and 
is in close interactions with the phytol chain ChID2 and the 
acyl chains of MGDG18, and to a lesser extent to SQDG4, 
DGDG6andMGDG7. 

The reduction/protonation mechanism of the secondary 
quinone has been a long-term challenge for researchers in 
photosynthesis. The Qb reduction occurs in two successive 
photoreactions including the formation of SQ intermediates. 
The first reaction is fast, occurring within 0.2-0.4 ms, and 
generates a SQ tightly bound into the Qb pocket, which can 
accept a second electron in a somewhat slower process - 0.6- 
0.8 ms [42]. Experimental data suggested that Qb' becomes 
protonated before the second electron transfer occurs, thus 
avoiding the formation of highly charged and unstable 
QaQb^ species [35, 36] and the protonation occurs more 
promptly if Qa is in its reduced state [10]. The resulting plas- 
toquinol is characterized by low affinity for the Qb binding 
niche and is promptly released into the membrane. 

In bRC the protonation process has been well estab- 
lished: the first H"^ necessary to neutralize the charge on the 
distal (in relation to the non-heme Fe position) keto-oxygen 
atom of the SQ is provided by the Ser L223 amino acid 
trough Asp L213, while the second one is released by the 
Glu L212 residue. In plant and cyanobacterial PSII RCs the 
details of this phenomenon are still under investigation. 
Similarly to bRC, it has been supposed that the first protona- 
tion is accomplished by Dl-Ser264 via Dl-His252, while 
different possibilities have been explored for the mechanism 
of the second protonation [35, 36, 43]. Recently, by large- 
scale quantum mechanical/molecular mechanical approach 
based on the 1.9 A resolution crystal structure, the Dl- 
His215 was suggested as the sought H'^ donor for the keto- 
oxygen of QbH proximal to the non-heme Fe [10]. The 
authors speculated that the resulting Dl-His215 anion, in 
turn, is protonated in a redox-process, which most probably 
includes bicarbonate, Dl-Tyr246 (Fig. 2, centre panel) and 
water as well as rearrangement of H-bonds network within 
the Qb binding niche. 

2.3. Additional PSII Plastoquinone Binding Sites 

A third plastoquinone molecule (Qc) associated to PSIIcc 

was identified in the 2.9 A resolution crystal structure of T. 
elongatus, and five of its isoprenoid units were successfully 
modelled into the PQ/PQH2 exchange channel I, between the 
TMH of PsbJ and Cyl.bssg a and jS-subunits [5]. In a subse- 
quent 3.6 A resolution crystallographic study of monomeric 
PSII core complex from T. elongatus, a trace of electron 
density was revealed in the same position but it was not pos- 
sible to assign it to a specific molecule [6]. No electron den- 
sity was detected in the putative Qc binding side in the 
most detailed 3D structure at 1.9 A resolution from 
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Table 2. Factors involved in PSII plastoquinone interactions within Qa, Qb and Qc binding niche, summarized from Hasagavra and 
Noguchi [23]. Amino acid residues providing attractive or repulsive effects with IAE'I>2.0 kcal/mol and lipids and cofac- 
tors providing attractive or repulsive effects with IAE'I>0.1 kcal/mol are presented. Amino acids, Upids and cofactors are 
ordered following the decrease of the interaction ener^ with the corresponding PQ molecule. 



Factors Involved in 
PQs Interactions 


PQ at Qa Binding Niche 


PQ at Qb Binding Niche 


PQ at Qc Binding Niche* 


Protein Matrix 


Head attractive interactions 


D2-Phe261,-Thr217, 
D2-Trp253,-His214, 
D2-Met246, -Ser262 


Dl-Phe255, -Phe265, 
Dl-Val219, -Leu218, 
Dl-His215 


None 


Tail attractive interactions 


D2-Gly200, -Gly203, 
Dl-Phe52,PsbL-Leu30, 
PsbT-PhelO, PsbL-Leu27, 
D2-Phe270,Dl-Ile77 


D2-Phe38, -Tyr42, 

Dl-Phe211, 
PsbF-Leu26, -Thr25, 
D2-Arg26, -Argl28 


PsbJ-Gly20, -Val21, 
PsbJ-Glyl8, -Alal7, 
PsbJ-Metl9, -Ile24, 
PsbF-Thr30,-Ile31, 
PsbE-Phe31 


Repulsive interactions 


Head: D2-Gly215, 

D2-Arg251 
Tail: D2-Leu279, 
D2-Val204, -Val274 


Head: Dl-Gly253, 
Dl-Asn266 
TaU: D2-Leu43, 

PsbF-Val23, -Ala22, 
Dl-Leu271,D2-Asp25 
PsbE-Glu7 


Head: PsbE 
Tail: PsbT 


Lipids and Cofactors'' 


Head atttactive interactions 




DGDG6 


MGDGIS, Chloz, MGDG7, Carnj, 
DGDG6, SQDG4 


Tail attractive interactions 




ChlD2, MGDG18, SQDG4, DGDG6, 
MGDG7 


CarD2, MGDG18, DGDG6 


Rcpulsi\ e inleraclions 




Head: MGDGIS. Chl„; 


Tail: MGDG7. Car,,;. 



"For Qc all PSII amino acid residues providing attractive effects with IAE'I>1.0 kcal/mol are presented. 

''The abbreviations for lipids and molecular cofactors are same as define in Section 1 and the List of Abbreviations. 



T. vulcanus [8], and this binding spot was neither occupied 
by PQ nor by PSII herbicide terbutryn in T. elongatus 
PSII/terbutryn crystal complex resolved at 3.2 A resolution 
[7]. The uncertainty of Qc presence in the different structures 
has been ascribed to differences in adopted purification and 
crystallization procedures, and photosynthetic sample pre- 
treatment, to the low binding affinity of Qc, or to the lack of 
a well-defined binding pocket [5, 7]. In fact, neither H-bonds 
nor 71-71 interactions were observed for Qc stabilization and, 
most probably, the molecule binding is a result of the possi- 
ble van der Waals interactions with the nearby hpid mole- 
cules, ChID2, CarD2 and phytyl tail of Qb [5, 23], Table 2, 
PQ Head attractive interactions and (Fig. 2). FMO calcula- 
tions based on the 2.9 A resolution structure pointed out that 
more than 90% of the Qc interaction energy with the protein 
matrix and with PSII integral lipids and cofactors is due to its 
isoprenoid chain (Table 1). In addition, the relative contribu- 
tion of the interactions with the surrounding protein matrix to 
Qc total binding energy is about 56%, that is much lower 
compared to that calculated for Qa and Qb [23] - about 
100% and 82%, respectively (Table 1). The high mobility of 
the Qc molecule has been attributed to its peculiar interaction 
with the environment, which is characterized by the head 
group rimming predominantly by a belt of hpids and cofac- 
tors, and the tail fencing to proteins (Table 2). The X-ray 



localization of a third PQ molecule in the 2.9 A PSII struc- 
ture supported earlier findings on the occurrence of addi- 
tional quinones (different from Qa and Qb), located in the 
PQ hydrophobic cavity close to the Cyl.bssg [44, 45]. The 
physiological role of the Qc is still under investigation. 
Guskov and coworkers [5] proposed three hypothetical 
mechanisms for PQ/PQH2 exchange, in which Qc is alternat- 
ing with Qb or is the PQ molecule waiting to enter into the 
QB-binging site, thus contributing to the rapid turnover of the 
PQH2 with a new PQ molecule from the membrane pool. 
Experimental evidences have suggested a strong involve- 
ment of Qc in a PSII photoprotective electron pathway medi- 
ated by Cytbssg [20-22]. It was hypothesized that Cyt./; 559 
could act as a plastoquinol oxidase (oxygen reductase) or a 
superoxide oxidase/reductase enzymes; these functions have 
been extensively documented in other reviews [15, 16, 36]. 
A distinctive feature of this heme-containing protein is its 
redox heterogeneity. In PSII membrane fragments 75% of 
Cyt.fejjs is in its high potential form, 16% is in its intermedi- 
ate potential form and 9% is in its low potential form; its 
redox properties and reactivity have been correlated to the 
occupancy of the Qc site [22]. A biphasic pattern of Cyt.bssg 
reduction by PQH2 has been observed in PSII membrane 
fragments containing predominantly oxidized Cyt.bssg, and 
experimental results suggested that a PQ/PQH2 molecule 
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Fig. (2). Top: Schematic view of PSII reaction centre showing the relative position of plastoquinones Q^, Qb and Qc as red sticks. Chloro- 
phylls, pheophytins and carotenoids are also shown (in green, blue and orange sticks, respectively). The OEC Mn4Ca cluster and the non- 
heme iron atom are represented by purple and red spheres, respectively. Dl and D2 backbone are shown as a blue and yellow ribbon, respec- 
tively. Bottom: Detailed view of the plastoquinones Qa (left panel), Qb (centre panel), and Qc (right panel) binding sites. Residues mentioned 
in the text are shown as sticks coloured by atom type. 



bound to the Qc site is responsible for the fast phase (with 
rate constant of 100 ms~') of the Cyl.bssg reduction [46]. The 
slow Cyt.bssg reduction phase (with rate constant of 2.4 
min~') was attributed to formation of redox equilibrium be- 
tween a set of Cyt.bssg redox forms and "free" PQ/PQHt 
from the membrane pool. The most likely mechanisms sug- 
gested by the authors for the Cyl.bssg fast reduction requires: 
/) a location preferably encouraging one-electron redox- 
reaction; //) a polar environment allowing accommodation of 
the PQ/SQ/PQH2 and occurrence of the so called L sub- 
stances in anionic form (such as ADRY reagents, dinoseb, 
high concentration of DCMU, tetraphenylboron); and ///) a 
fine-tuned Cyt.iijjs/quinone interaction requiring well- 
defined quinone-protein interactions. However, the structure 
of the Qc binding site described in the 2.9 A resolution crys- 
tallographic study does not fulfil any of these requirements. 
The discrepancy between the structurally identified and ex- 
perimentally suggested additional PSII PQ binding sites was 
found also in previous studies [20-22] . In order to distinguish 
between the Qc site and the presumed polar quinone site of 
PSII that binds PQ/SQ/PQH2/L, the authors proposed to 
assign the latter as Qd [46]. These experiments led to hy- 
pothesize the existence of two different PQ binding sites, Qc 
and Qd, which host two different PQ molecules or one bound 
quinone that can occupy two different positions. 



3. DYNAMICS STUDIES IN PHOTOSYNTHESIS RE- 
SEARCH 

3.1. Neutron Scattering Techniques 

The determination of protein structure and the analysis of 
conformational changes are the key necessary to understand 
the biological function and the activity of a protein. The 
structural characterization of proteins, however, constitutes 
only part of the puzzle because their function requires a spe- 
cific dynamic. The protein motions arise on a wide range of 
timescales, from femtoseconds for bonds vibrations to milli- 
seconds or seconds, for more complex processes that require 
a structural rearrangement. The dynamics of proteins con- 
tributes to the thermodynamic stability of their functional 
states and plays an important role in various processes, such 
as allosteric processes, where the correlations between struc- 
tural fluctuations are able to transmit information between 
distant sites of the protein. 

The study of dynamic processes at the microscopic level 
has been developed using a wide range of high resolution 
spectroscopies (Mossbauer spectroscopy, nuclear magnetic 
resonance (NMR), time-resolved fluorescence). We should 
add elastic, quasi-elastic and inelastic neutron scattering 
(ENS, QENS and INS, respectively) technique, which offers 
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the advantage of covering length/distance (0.5-10 A) and 
energy (0-200 meV) scales that correspond to thermal fluc- 
tuations and provide information on the geometry of the 
movements [47, 48]. The neutrons as a dynamics probe are 
highly sensitive to movements of hydrogen atoms, which are 
uniformly distributed within biological structures and repre- 
sents nearly half of the atoms in them. On temporal and spa- 
tial scales accessible with neutron scattering, the hydrogen 
atoms trace the movements of the chemical groups, which 
they are linked to, such as amino acids side chains. In addi- 
tion, it is also possible to mask part of the signal using spe- 
cific deuteration, since deuterium has a weak scattering sig- 
nal with respect to that of hydrogen [48]. The neutron scat- 
tering experiments, with their peculiar chance to exploit the 
isotopic substitution technique, permit to integrate and give 
complement to data obtained by numerical simulation and 
through specific molecular biology techniques. Moreover, 
they provide unique information useful for the characteriza- 
tion of the dynamical heterogeneity of biological macro- 
molecules, from globular protein up to membrane com- 
plexes, in the temporal range 10-'^ - 10-' s. 

The neutron scattering is an invaluable technique to 
tackle a plethora of important open questions in Biophysics 
such as the definition of correlations between structure, dy- 
namics and functionality, the correlations between the inter- 
action of water molecules with proteins and its role in the 
folding process and enzymatic activity, the observation of 
how the hydrogen bond network between solvent and bio- 
logical macromolecules is affected and modified by selective 
genetic mutation, and hence the local flexibility and biologi- 
cal efficiency of the molecules [49-51]. 

3.2. Neutron Scattering in Photosynthetic Research 

Photosynthetic research exploring the relationship be- 
tween the function and the dynamics is a new field of inter- 
est. However, neutron scattering has been demonstrated to be 
a suitable tool to investigate protein dynamics in both com- 
plex multi-subunit photosynthetic systems and whole cells 
[52, 53]. 

A pioneer quasi-elastic neutron scattering study on the 
purple bacterium Rhodobacter sphaeroides RC proteins has 
demonstrated the effect of genetic point mutations for the 
overall protein dynamics [52]. The RC dynamics over a tem- 
perature range of 280 K was probed in a wild type and two 
non-functional mutants, carrying alanine residues on the 
place of Glu212 and Asp213 in the L-subunit. The Glu212 
residue is highly conserved in bRC and is involved in the 
secondary ubiquinol protonation (see section 2.2). The crys- 
tallographic structure of the mutant proteins showed a local 
configuration of the active site more open than that of the 
wild type. One of the parameters indicative for the internal 
protein motions, the mean square displacement (<u^>), cal- 
culated for the two mutants resulted to be higher than for the 
wild type protein suggesting a higher flexibility [52] . These 
results correlated the enhanced flexibihty of mutated RC 
proteins with a lack of functionality, thus suggesting the re- 
quirement of a rigid core for the accomplishment of protein 
function in a multi-protein complex, such as RC involved in 
highly tuned photosynthetic reactions. This ground-breaking 
experiment opened the way to the urgency to understand if 
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the core rigidity is essential for the charge transfer process, 
and if this property is shared by all the photosynthetic sys- 
tems, and how this information can be applied to the design 
of bio-sensors and organic semiconductors. 

An original green microalga species from the Chlorococ- 
cal order, which was never studied for its high tolerance to 
radiations, has been isolated in the Institute Laue Langevin 
research reactor storage pool. The microorganism able to live 
and reproduce in extreme environments was investigated by 
Farhi and coworkers combining neutron scattering experi- 
ments of internal dynamics with results of radiation effects 
on the microalgae physiology, metabolism and genomics 
[53]. The microalgae were stressed with radiation doses up 
to 20 kGy (2 Mrad), and studied by NMR, looking for modi- 
fication in the metabolism, and by ENS, looking for both 
dynamical flexibility and structural macromolecular changes 
in the cells. While the sugar metabolism was slightly reduced 
and only at high radiation dose, it was probed that the effect 
of a high irradiation on micro-algae does not lead to protein 
denaturation, but rather to reduction of macromolecules 
flexibility - as showed by the reduction of mean square dis- 
placement (<u^>) in (Fig. 2) of reference [53] . 

Pieper and coworkers have performed exhaustive neutron 
scattering studies on protein dynamics of PSII membrane 
fragments and antenna complex LHCII as a function of the 
temperature and sample hydration level [54]. They showed 
that PSII membrane fragments undergo a hydration depend- 
ent "dynamical transition" at about 240 K, underling the cru- 
cial role of hydration water for PSII flexibility. Above this 
temperature, PSII membrane fragments exhibit structural 
flexibility due to the presence of fast (picosecond) internal 
protein motions, while there is only harmonic vibrational 
dynamics at lower temperatures. It was then speculated that 
diffusive protein dynamics are indispensable for enabling Qa 
reoxidation by Qb at temperatures above 240 K, which ex- 
plains the strong dependence of this electron transfer step on 
temperature and hydration level of the sample. 

A newly developed technique of time-resolved QENS of- 
fers the advantage of simultaneous monitoring of protein 
internal dynamics and functional activity in specific sample 
functional states. This method was firstly applied for dynam- 
ics-function correlation in bacteriorhodopsin, a membrane 
embedded protein that functions as a light-driven proton 
pump in Halobacterium salinarum, and lately for the more 
complex PSII membrane fragments with inhibited Qa to Qb 
electron transfer [54, 55] . At physiological temperatures the 
calculated mean-square displacement of the hydrogen atoms, 
characterizing the flexibihty of the system, revealed a sig- 
nificantly increased flexibility between 310 and 320 K [56]. 
The authors claimed that this transition is evidently corre- 
lated with the detachment of OEC from the membranes, as 
suggested by its absence in samples in which OEC has been 
removed prior to the QENS experiments. The transition is 
not accompanied by significant changes in the molecular 
organization of the pigment protein complexes, but the role 
of monomerization of PSII dimers cannot be ruled out. 

Lately, preliminary results on intracellular water collec- 
tive dynamics investigation on Chlamydomonas cells, carry- 
ing both native and mutated Dl RC proteins, have been also 
addressed. The authors proposed a comparison between fuUy 
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hydrated and partially de-hydrated deuterated whole cells. It 
was observed a distinct sound propagation speed between the 
two hydration levels suggesting a more rigid structure of 
hydration water than intracellular water [57]. 

4. COMPUTATION ANALYSES FOR MODIFICA- 
TION OF Qb binding pocket amino acid 

COMPOSITION: FUNDAMENTAL RESEARCH AND 
APPLICATIONS 

In the last decades, several mutational studies on Dl pro- 
tein were conducted to determine the effects induced by 
aminoacidic substitutions located in the stromal loop of the 
D 1 protein. Mutations resulted in an impairment of the elec- 
tron transfer between Qa and Qb leading to the identification 
of residues often involved in modulation of PSII photoinhibi- 
tion susceptibility and/or herbicide specificity [58-62]. Com- 
putational studies on PSII Dl protein aimed at rational modi- 
fication of the Qb binding pocket have been scarce partly 
because a high resolution three-dimensional structure of PSII 
has been available only relatively recently and partly due to 
the difficulties in generating site-specific mutants of this 
protein. 

The first structure of PSII from T. elongatus has been 
published only in 2004 and at a resolution of only 3.5 A 
[PDB code 1S5L, 3]. The same year a 3.2 A resolution struc- 
ture has been published independently [PDB code 1W5C, 
63] and refined the following year at 3.0 A resolution [PDB 
code 2AXT, 4]. However, this last structure had still an in- 
sufficient resolution to uncover all the atomic details of this 
macromolecular proteins-cofactor complex such as the pre- 
cise architecture of the OEC. Nonetheless, the Loll et al. 
crystal structure was sufficiently detailed to allow the first 
computational studies on the Qb binding pocket as well as 
homology modelling of PSIIs from other organisms. The 
first of such studies was probably that of Loll and coworkers 
[64] who modelled the three different copies of T. elongatus 
Dl protein (PsbAl-A3) within the 3.0 A resolution crystal 
structure of PSII, revealing that most of the Dl amino acid 
variants were located in the direct vicinity of redox-active 
cofactors of the electron transfer chain. However, proper 
computational studies aimed at selecting amino acid variants 
within the Qb binding pocket able to modify the recognition 
properties of the pocket itself appeared only the following 
year [65]. In this study, homology modelling of the Chlamy- 
domonas reinhardtii Dl protein was carried out based on the 
T. elongatus crystal structure, followed by in silico mutage- 
nesis and energy minimization analyses to identify functional 
amino acids in the Dl protein interacting with herbicides. 
Following the results of theoretical calculations, three site 
directed mutants were produced by site-directed mutagenesis 
and characterized by fluorescence analysis. This study dem- 
onstrated that Dl mutants S268C and S264K could be used 
in the development of biomediators suitable for the selective 
detection of triazine and urea classes of herbicides, respec- 
tively. The study was further refined the same year undertak- 
ing an in silico study of the C. reinhardtii D1-D2 proteins 
with the aim of designing mutants with increased affinity for 
atrazine [66]. The final outcome of this work was the design 
and selection of Dl mutants hypothetically able to increase 
atrazine binding affinity by orders of magnitude, represent- 



ing a useful theoretical guide for the development of high 
affinity atrazine biosensors. 

With the increase in computing power, it has become fea- 
sible to adopt more sophisticated structural bioinformatics 
approaches, such as molecular dynamics (MD) simulations, 
for the study of the recognition properties of the Qb binding 
pocket toward natural and exogenous ligands. The first all 
atoms MD simulations of the complete T. elongatus PSII 
structure (including all cofactors) embedded in a lipid mem- 
brane was published in 2011 [67]. The aim of this work was 
that of uncovering the fine details of atrazine molecular in- 
teractions within the Qb binding pocket using a full atomistic 
(and thus realistic) description of PSIIcc. Atrazine was first 
docked in the vicinity of the Qb binding pocket and then let 
free to diffuse during a 10 ns MD simulations trajectory. 
During the simulations atrazine bound deeper into the Qb 
pocket establishing hydrophobic interactions with PSII Dl- 
Phe255 phenyl ring and Dl-Met214 side chain, and a single 
strong hydrogen bond with Dl-His215 (Fig. 3). On the other 
hand, no binding partner appeared to stabilize other groups 
of the atrazine molecule such as the chloride atom (bound in 
an energetically unfavourable position in the vicinity of the 
aliphatic residue Dl-Val219) and the N4 and N5 protons for 
which no hydrogen bonding partner was observed. The con- 
clusions of this study were that rational design of site- 
directed mutants was likely to increase atrazine affinity for 
the Qb binding pocket. In particular, it was suggested that, 
among others, mutation of Dl-Phe265 residue with a polar 
hydrogen bond acceptor residue could stabilize the atrazine 
N5 proton and atrazine binding within the pocket. 




Fig. (3). Schematic view of the interactions of atrazine within the 
PSU Qb binding pocket predicted through a 10 ns MD simulations 
run [66]. For clarity, only the Dl backbone (orange ribbon) and 
residues mentioned in the text are shown. The dashed green line 
highlights the hydrogen bond between atrazine and His215. 
Atrazine carbon atoms are in green, nitrogen atoms are in blue, 
sulphur atom is in yellow and the chloride atom is in magenta. 

CONCLUSIONS AND PERSPECTIVES 

The advances in the crystallographic characterization of 
PSIIcc have broaden the understanding of PSII function and 
organization allowing to identify the atomic interactions tak- 
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ing place between the RC protein matrix, electron tremsport 
cofactors and thylakoid membrane lipids. Unravelling finer 
molecular details and mechanisms of the PSII structure, 
function and dynamics is mandatory to support novel bio- 
technological applications exploiting its unique properties of 
charge separation, electron transfer and water splitting. Spe- 
cifically, it is possible to exploit PSII enriched particles or 
their active subcomponent to develop receptor elements suit- 
able to be exploited in photosynthesis-based biosensors for 
the detection of environmental pollutants. The latter, in fact, 
can interact with the PSII donor side interrupting the light- 
induced electron transfer, generating a measurable opto- 
electronically detectable signal. In this context, insights into 
atomic interactions and dynamics of polypeptides and co- 
factors provided by both molecular dynamics or docking 
studies will help to design novel bio-recognition elements 
having improved selectivity and sensitivity towards different 
classes of compounds. Studies exploring the relationship 
between the function and the dynamics of the photosynthetic 
complex systems are still in their nascent stages, along with 
MD simulations. However, these techniques could offer a 
real possibility of investigating the dynamics associated with 
a molecule's biological function and complement the pro- 
gress of this technology. Besides, the extraordinary proper- 
ties of water- splitting could be used to design bio-inspired 
device for bio-energy production, while the charge separa- 
tion features could be exploited to design molecular photo- 
voltaics. 
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LIST OF ABBREVIATIONS 



ADRY - Acceleration of the deactivation reactions 

of system Y 

bRC = Bacterial reaction centre 

Car = Carotenoids 

Chi = Chlorophyll 

CP43 = PSIIcc chlorophyll apoprotein encoded by 

psbQ. gene 

CP47 - PSllcc chlorophyll apoprotein encoded by 

psW> gene 

Cyt.b559 = Cytochrome h^^g 



DGDG = Digalactosyldiacylglycerol 

ENS = Elastic Neutron Scattering 

FMO = Fragment Molecular Orbital 

INS = Inelastic neutron scattering 

L substances = Chemical compounds that induce negative 
shift in the midpoint potential of high po- 
tential form of Cyt./7jj9 

LHCll = Light Harvesting Complex II 

MD = Molecular Dynamics 

MGDG = Monogalactosyldiacylglycerol 

NMR = Nuclear Magnetic Resonance 

OEC - Oxygen Evolving Complex 

PDB - Protein Data Base 

PG = Phosphatidylglycerol 

RC - Reaction Centre 

Pheo = Pheophytin 

PQH2 = Plastoquinol 

PQs = Plastoquinones 

PsbE and 
PsbF PSIIcc 

proteins = Corresponding to a- and j5-subunit of cyto- 
chrome ^559 

PSI - PhotosystemI 

PSII = Photosystem II 

PSIIcc = PSII core complex 

Qa - PSII primary quinone electron acceptor 

Qb = PSII secondary quinone electron acceptor 

Qc and Qd = PQ molecules supposed to be bound to 
PSIIcc 

QENS = Quasi-elastic neutron scattering 

SQ - Semiquinones 

SQDG = Sulfoquinovosyl-diacylglycerol 

terbutryn = 2-ferf-Butylamino-4-ethylamino-6- 
methy Ithio- 1 ,3,5-triazin 

TMH = Trans-membrane helices. 
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